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The spin-Peierls instability describes a structural transition of a crystal due to strong magnetic interactions. 
Here we demonstrate that cold Coulomb crystals of trapped ions provide an experimental testbed in which 
to study this complex many-body problem and to access extreme regimes where the instability is triggered 
by quantum fluctuations alone. We present a consistent analysis based on different analytical and numerical 
methods, and provide a detailed discussion of its feasibility on the basis of ion-trap experiments. Moreover, we 
identify regimes where this quantum simulation may exceed the power of classical computers. 
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The beauty of low-dimensional quantum many-body sys- 
tems (QMBS) lies in the complexity born of the combination 
of interactions, disorder, and quantum fluctuations. However, 
these ingredients also conspire to render perturbative tech- 
niques inefficient, posing thus a fundamental challenge that 
has inspired the development of a variety of analytical [1] 
and numerical [2] tools. Moreover, the synthesis of low- 
dimensional materials has upgraded these challenges from 
a theoretical endeavor into a discipline that underlies some 
of the most exciting recent discoveries in condensed-matter 
physics, such as the fractional quantum Hall effect. The recent 
progress in the field of atomic, molecular, and optical (AMO) 
physics presents a promising alternative to these solid-state 
realizations of low-dimensional QMBS. This field, which was 
originally devoted to the study of light-matter interactions at 
the scale of a single or few atoms, is progressively focusing on 
the many-body regime in platforms such as neutral atoms in 
optical lattices [3], cold Coulomb crystals of trapped ions [4], 
or coupled cavity arrays [5]. The possibility of experimen- 
tally designing the microscopic Hamiltonians in order to tar- 
get a variety of complicated many-body models introduces a 
novel approach to explore QMBS in a controlled fashion, the 
so-called quantum simulations (QSs) [6]. Some remarkable 
QSs in AMO platforms are the optical-lattice realization of 
Mott [7] and Anderson [8] insulators, and the recent efforts 
towards a quantum-Hall insulator [9]. More specific to the 
subject of this manuscript is the QS of quantum magnetism in 
trapped ions [10], and more recently in optical lattices [11]. 

In this Letter, we explore the capabilities of AMO setups 
for the QS of interaction-mediated instabilities in QMBS. The 
standard playground for these phenomena is the so-called 
metal-insulator transition [12], which has been investigated 
for a variety of transition-metal compounds in the field of 
strongly-correlated electrons. A paradigmatic case is the one- 
dimensional metal, where either the electron-electron interac- 
tions destabilize the metal towards a superconducting state, or 
the electron-phonon coupling leads to a charge-density-wave 
condensate [35]. The latter instability is a consequence of the 
so-called Peierls transition [14], where the electron-phonon 
interactions induce a periodic distortion of the ionic lattice, 
and open an energy gap in the conduction band of the metal. 
By virtue of the Jordan-Wigner transformation [10], this phe- 
nomenon finds a magnetic counterpart: the spin-Peierls tran- 
sition [8], whereby a spin-phonon-coupled antiferromagnet 
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Figure 1 . Scheme of the spin-Peierls instability: (a) In the para- 
magnetic phase |P) = | tt " ■ " t) a H the spins are parallel to a trans- 
verse field g> g c that points along the z-axis, and the Coulomb crys- 
tals corresponds to an ion string, (b) The antiferromagnetic phase 
g < gc corresponds to the two Neel-ordered groundstates |AF) G 

{| H 1 — ), | — | h)}, where the spins are antiparallel in 

the x-basis |±) = (|t)±|i)) / y/2. This order-disorder quantum phase 
transition occurs with the linear-to-zigzag structural phase transition, 
(c) Arrangement of the laser wavevector kL lying within the xy-plane. 



becomes unstable with respect to a dimerization of the lattice. 
This creates an alternating pattern of weak and strong spin 
interactions, which in turn opens an energy gap in the spec- 
trum of collective excitations. We note that this instability has 
turned out to be important for different compounds, such as 
organic molecular crystals and transition-metal oxides [17]. 

From a theoretical perspective, the complete understand- 
ing of such a complex many-body system, treating the dy- 
namics of the spins and phonons on the same footing, is still 
considered to be an open problem [17, 18]. From an exper- 
imental point of view, the spin-Peierls instabilities observed 
so far [17, 19] take place at finite temperatures, and are lim- 
ited to the so-called Heisenberg model. Hence, the possibility 
of realizing a spin-Peierls transition only driven by quantum 
fluctuations, and possibly exploring different types of mag- 
netic interactions, remains as an experimental challenge. 

We hereby present a theoretical proposal for a trapped- 
ion QS to tackle both problems. In particular, by building 
on the recent experiments [10] on the quantum Ising model 
(QIM) [1 1], we describe how to tailor a spin-Peierls instabil- 
ity. We show that (i) the disordered paramagnet in a linear 
ion chain changes into an ordered antiferromagnet in a zigzag 
crystal [Fig. l(a)-(b)], and (ii) the spin-Peierls transition can 
be driven only by the quantum fluctuations introduced by the 
transverse field of the QIM. Let us remark that, in compar- 
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ison to neutral atoms in optical lattices or coupled arrays of 
cavities, trapped ions seem to be the best candidates to realize 
the spin-Peierls quantum simulator. One of the main reasons 
is that the underlying lattice is not externally fixed, but rather 
results from the self-assembling dynamics of the ions. 

The system- The advent of experimental techniques for the 
confinement, cooling, and coherent manipulation of atomic 
ions has recently been exploited for QS purposes [4], where 
the controlled increase of the number of trapped ions yields 
a genuine bottom-up approach to the many-body regime. We 
consider a Coulomb gas formed by an ensemble of N trapped 
ions of mass m, and charge e, which are described by the 
Hamiltonian 

(1) 

where {cQa}a=x,y,z are the effective trapping frequencies of a 
linear Paul trap. Here, COo is the energy difference between 
two electronic groundstates of the atomic structure | t)*, | i)i, 
where ft = 1, and of = | t/)(t/| _ I U)(U\- Tms Hamiltonian 
must be complemented by the laser-ion interaction responsi- 
ble for coupling the electronic and the motional degrees of 
freedom. We consider a pair of laser beams with frequencies 
{c0/}/ = i 5 2, wavevectors {k/}/ = i2, and phases {0/}/=i,2> which 
are tuned close to the atomic transition. In the dipolar approx- 
imation, the laser-ion Hamiltonian becomes 

H L = £(£1/C7+ + flfar) cos(k/ • r,- - co t t + 0/), (2) 
i,i 

where £1/ stands for the Rabi frequency of the transition, and 
we have introduced the spin raising and lowering operators 
a t = I ti) (U\ — ( G i~y • T° proceed further, we need to make 
some assumptions about the dynamics of this atomic plasma. 

As evidenced in early experiments [22], a laser-cooled en- 
semble of ions self-assembles in a Coulomb crystal, which 
undergoes a series of structural phase transitions (SPTs) as 
the trapping conditions are modified. In particular, when 
C0y ^> co x ,co z , and the ratio k x = (co z /co x ) 2 is tuned across a 
critical value Kq [5], the geometry of the crystal changes from 
a linear string to a zigzag ladder. Note that this SPT displays 
a rich phenomenology that has recently revived the interest 
in the subject [1, 4, 6, 25, 26, 29]. Here, we focus on the 
linear regime close to the critical point K x < K^, where the 
vibrations of the ions along each of the confining axes are 
decoupled. We consider that the laser wavevectors in (2) lie 
within the xy-plane [Fig. 1(c)], k/ = k[ x e x + ki y e y , such that 
their frequencies are tuned close to the resonance of the vi- 
brational sidebands (0/ w COo =•= (Oy> Since these correspond to 
the strongly-confining y-axis, C0y ^> (0 X , co z , the coupling of the 
laser beams to the x and z vibrational modes becomes far off- 
resonant and can be neglected. Let us remark that this argu- 
ment has one possible exception, there might be a vibrational 
soft mode where the phonons condense at the SPT. As identi- 
fied in [1], this corresponds precisely to the zigzag mode along 
the x-axis. This affects the laser-ion coupling (2) regardless of 
the soft-mode frequency (see Supplementary Material). 

The above considerations allow us to extract the relevant 
part of the Hamiltonian (1) after introducing r ; = l z (^e z + 



qix^x + qiy£y + Q iz^z)-> wnere ^ are the equilibrium positions 
in units of l z = (e 2 /mco 2 ) 1 / 3 , and qi a are the small displace- 
ments along the corresponding axes. Following [6], the dis- 
placements along the x-axis have been adapted to the afore- 
mentioned zigzag mode qi x = (— l)'<5g iX , such that Sqt x is a 
smooth function that allows a gradient expansion. The Hamil- 
tonian then becomes Hq = \ COotf + H x +H y , where 

(3) 

and we have introduced the gradient djfj = fi — fj. In these 
expressions, the coupling energies for the vibrations are 

ri=mco^l 2 z (1 - 1*^(3)) , k? = «fl& 2 (^£(5)) , 

along the x-axis, expressed in terms of the inhomogeneous 
function Q(n) = Ez#[(-l) f - (-l) 1 }"- 1 |z° - ^\~ n with n £ 
Z, and the Kronecker delta <5/ m . In the limit of tight confine- 
ment along the y-axis, K y = (co z /cOy) 2 <C 1, we find 

i o m(D 2 l 2 K y 

r? mco 2 l 2 K?. y (5) 

Accordingly, H y corresponds to a set of dipolarly-coupled 
harmonic oscillators, whereas H x describes a set of nearest- 
neighbor-coupled anharmonic oscillators. 

The coupled harmonic oscillators (3) are diagonalized 
yielding a set of collective modes with frequencies co n , whose 
excitations, created and annihilated by a^a n , shall be referred 
to as the hard phonons. This yields the quadratic Hamilto- 
nian H y — Y,n dO n a\a n . By setting the laser frequencies to the 
red and blue vibrational sidebands of the atomic transition, 
C0\ ~ COo co n and o>2 « COo + co m we can express the laser-ion 
interaction (2) as follows (see Supplementary Material) 

Hl = £ (& r J e ^c+a n + ^lt^a+al + H.c.) , (6) 

in 

where we have introduced the sideband coupling strengths 
&\ n = \a X T\\nJtin^ x = \£lir[ 2nJ(in^ \ the collective 
Lamb-Dicke parameters r\i n — ki y /^/2mco n <C 1, the normal- 
mode vibrational amplitudes <Mi n , and 6 r = k\ x l ZJ % = kixh- 

In contrast, the anharmonic oscillators (3) correspond to an 
inhomogeneous version of the (f) 4 model on a lattice, namely, 
an interacting scalar field theory that cannot be exactly diag- 
onalized. This model yields a SPT that can be understood as 
follows. The regime rf > 0, u* > corresponding to trapping- 
frequency ratios fulfilling K x < = 2/Q(3), yields the lin- 
ear ion configuration, which respects the Z2 symmetry of the 
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model. Conversely, when rf < 0, u* > for K x > /c c ; , the 
ions self-organize in the the zigzag ladder corresponding to 
the broken- symmetry phase, whereby (Sqi x ) ^ signals the 
condensation of the soft phonons in the zigzag mode. We note 
that the (j) 4 model fulfills rf ^ r^, which leads to an inhomo- 
geneous SPT setting at the center of the trap [26]. As out- 
lined previously, when the soft phonons condense (Sqix) 0> 
there is a non-trivial effect in the laser-ion Hamiltonian (6) 
that must be considered. We show below that, in this case, 
the hard phonons mediate a spin-spin interaction, whereas the 
soft condensed phonons are responsible for a dimerization of 
the coupling strengths. This turns out to be the key ingredi- 
ent for a zero-temperature spin-Peierls transition. Let us also 
emphasize that this model, which is a cornerstone in the mi- 
croscopic description of SPTs [30], has not been combined 
with a spin-Peierls distortion to the best of our knowledge. 

Dimerized quantum spin model - The spin-phonon model 
in (3) and (6) yields an extremely complex QMBS. We ana- 
lyze the onset of a spin-Peierls quantum phase transition by 
performing a series of simplifications. First, we neglect the 
time-dependence of the zigzag distortion (3). This adiabatic 
approximation, which is standard in the treatment of spin- 
Peierls phenomena [8], is justified if the zigzag mode is much 
slower than the effective spin dynamics, which is valid close 
to the critical point. Hence, we treat the SPT classically by set- 
ting qi x = (—l) l (Sqi x ) self-consistently. Second, we consider 
a homogeneous zigzag distortion, which amounts to neglect- 
ing the nearest-neighbor couplings in Eq. (3). Third, when 
the coupling of the spins to the hard phonons (6) is weak, 
they can be integrated out yielding an effective quantum spin 
model. In the linear string, this leads to a dipolar version of 
the celebrated QIM [31], whereas the frustrated J1-J2 QIM 
arises in the zigzag configuration [32]. In this work, we show 
that in the vicinity of the critical point K x « Kcj, the quan- 
tum spin model corresponds to a dipolar QIM with additional 
spin-spin couplings whose sign alternates periodically when 
the soft phonons condense. 

In analogy to the S0rensen-M0lmer gates [33], we consider 
that the red- and blue- sideband terms (6) have opposite de- 
tunings S m — —S n \) =: Sn, where 8 m = (0\ — (cdq — co n ), and 
S n b — °>2 — (®o + G*n)- Besides, their Rabi frequencies fulfill 
£21*^ = ^kjy, and attain values such that \^ n \ = <C S n . 
In this limit, the sidebands (6) create a virtual hard phonon 
which is then reabsorbed by a distant ion, leading thus to 
the aforementioned effective spin-spin interaction. The above 
constraints are responsible for the destructive interference of 
the processes where a phonon is created and then reabsorbed 
by the same ion, a crucial property that underlies the availabil- 
ity of an effective spin Hamiltonian that is decoupled from the 
hard phonons. Finally, by considering that the pair of laser 
beams are counter-propagating ki = — k2 =: k, which implies 
that the parameters are 6 r = — % =: 0, it is possible to obtain 
the following spin Hamiltonian (see Supplementary Material) 



where the coupling strengths are the following 

Jfj = Jij ( cos ( 6 (q ix - q jx )) + cos ( 6 (q ix + q jx ) + 0_ ) ) , 
J IJ = Jij{cos(6(q ix - q jx )) - cos (6(q ix + q jx ) + 0_) ) , 
J ij = J ij (sin(0 (q ix -qjx))- sin(0 (q ix + q jx ) + 0-)), 
J ij = ~ J ij ( sin ( (to - 1jx) ) + sin( (q ix + q jx ) + 0- ) ) , 

and we have introduced the relative phase between the lasers 
(f)- = 0i — 02- Here, the spin-spin coupling strengths are 

where we made a expansion for ^ « 1, and introduced the 
bare detuning 8 y = C0\ — (coo — C0y) = — o>2 + + ®y), the 
bare Lamb-Dicke parameter r} y = k y / ^/2m(Oy, and the com- 
mon Rabi frequency := ^1 = ^2- Hence, the spin-phonon 
interaction leads to a generalization of the famous XY quan- 
tum spin model [7]. In order to obtain the promised dimerized 
spin model, one has to phase-lock the laser beams 0_ = 0, 
and consider the critical region where the zigzag distortion is 
small enough 0(8 qi) <C 1. Then, we can Taylor expand and 
obtain an antiferromagnetic Ising interaction characterized by 
Jff = 2Jij, and jJJ = 0, but also 

J% = 2J ij (-iy+ 1 0(8q jx )=J%, (10) 

which give rise to the quantum dimerization (i.e. a magnetic 
interaction that does not commute with the Ising coupling, and 
alternates between a ferromagnetic-antiferromagnetic sign 
jfjiJj*oz (— l) i+1 ). Besides, we also consider a transverse 
field h that can be obtained from a microwave that is far off- 
resonant with respect to the atomic transition. Altogether, the 
spin Hamiltonian becomes // e ff = Hising +#dimer, where 

"dimer-L J ||0_g0|3 ^ ' 

and we have introduced ^ = 0(5^) <C 1. Let us remark that 
the couplings of the dimerization Hamiltonian depend on the 
condensation of the soft phonons (Sqix), which is in turn de- 
scribed by the (j) 4 theory. Below, we show how this model 
leads to the desired Spin-Peierls instability. 

Spin-Peierls quantum phase transition - To demonstrate 
that the introduced scheme yields a QS of the spin-Peierls in- 
stability, we simplify the model by neglecting its long-range 
interactions and inhomogeneities. The spin model can be 
solved by means of a Jordan-Wigner transformation after set- 
ting £ = V/, which is later used as input to the (j) 4 model 
self-consistently. In the limit £ <C 1, the groundstate energy 
fulfills 

£ g (§)«£ g (0)-^9te)« 2 <^g(0), (12) 
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Figure 2. Spin-Peierls transition: (a) Scaling of the ground-state 
energy with the lattice dimerization. (b) Displacement of the critical 
point calculated from the divergence of the magnetic susceptibility. 



where / > is the nearest-neighbor antiferromagnetic cou- 
pling, g = h/J, and we have introduced a monotonically- 
decreasing positive-definite function <p(g) that depends on the 
complete elliptic integrals (see Supplementary Material). We 
have compared this expression to numerical DMRG [14] cal- 
culations (see Fig. 2(a) and Supplementary Material), which 
support our claim for small dimerizations. The above lower- 
ing of the groundstate energy pinpoints the instability towards 
the lattice distortion. Besides, the spectrum of magnetic exci- 
tations displays the following energy gap A oc |g — ^\ -f 4^ 2 | . 
With respect to the paramagnetic-to-antiferromagnetic quan- 
tum phase transition of the standard QIM at g c = 1 [11], the 
dimerization breaks the self-duality of the model and allows 
the critical point to flow towards higher values 



(13) 



In Fig. 2(b), we obtain the critical point from the divergence 
of the magnetic susceptibility Xm ^ Xm = —d 2 E g /dg 2 , which 
clearly shows that the critical point g c = 1 flows towards 
higher values g c (%) as me dimerization increases. 

Therefore, the disordered paramagnet close to criticality 
g c < g < g c will be unstable towards the ordered antiferro- 
magnet if the lowering of the energy (12) compensates the 
structural change. For self-consistency, we incorporate this 
energy change in the (f) 4 model. By noticing that it fulfills 
E g (^) — E g (0) oc £• Sqf x , it becomes clear how to modify the 
parameters of the scalar field theory (4). In analogy to the 
magnetic quantum phase transition, the critical value of the 
SPT is also displaced, but this time towards a smaller value 



C«(3) , 27 2 <p(g) 



ma, 2 /, 2 



K 



(14) 



Therefore, the linear ion string close to criticality k Ci i < K x < 
Kcj is unstable towards the symmetry-broken zigzag phase. 

We have thus proved our claim (i) that the paramagnetic 
phase in the linear ion string will be unstable towards the anti- 
ferromagnetic zigzag ladder [Fig. 1]. Moreover, by fixing the 
ratio of the trapping frequencies in the linear regime, K x < K C j, 
we can drive both the structural and the magnetic phase transi- 
tions by only modifying the transverse magnetic field g across 
g c . The necessary condition is that the trapping frequencies 



are tuned according to the following expression 

270 2 <p(fc) 



mCO z l z 



71 fa 1 - K*'( 3 )) ' 



(15) 



Hence, the zigzag antiferromagnet g < g c is driven onto a lin- 
ear paramagnet by increasing the quantum fluctuations g > g c 
(and vice versa). This supports our claim (ii) that the spin- 
Peierls transition can be driven by quantum fluctuations alone. 

Experimental considerations - We focus on 25 Mg + and se- 
lect two hyperfine levels for the spin states | ti) = \F = 2, mp = 
2), | li) = 1 3, 3), such that the resonance frequency in (1) is 
(Oq/2k = 1.8 GHz. We consider a N = 30 ion register with 
trapping frequencies (O z /2n « 300 kHz, (O x /2n « 7 MHz, 
and COy/271 = 10 MHz. The phase-locked laser beams lead- 
ing to (2) are blue-detuned Sy/2n ^1 MHz, such that the 
two-photon Rabi frequencies are £1^/271 w 1 MHz, and the 
Lamb-Dicke parameter r\ y « 0.2. With these values, the re- 
quired constraints detailed in the Supplementary Material are 
fulfilled, and we obtain a nearest-neighbor spin coupling with 
the typical strength / = 2/ e ff/|z? — 1 3 « 1 kHz observed in 
experiments [10]. By considering these parameters, the condi- 
tion (15) imposes the following constraint over the anisotropy 
(Kqj — K x )/Kcj ~ 10 _4 7] 2 , which requires to be sufficiently 
close to the structural phase transition. In practice, the soft ra- 
dial trapping frequency must be controlled with an accuracy of 
Acd x ~ 10 -6 co x ^1-10 Hz, which coincides with the the preci- 
sion required to observe quantum effects in the SPT [6] . Pro- 
vided that this precision is achieved in the experiments, one 
could optically pump the linear ion register to | \jf(0) ) = <£> | ti) > 
and then study its adiabatic evolution towards the AF phase as 
the transverse field g(t) is decreased. The corresponding AF 
order can be measured by fluorescence techniques, whereas 
the structural phase transition could be directly observed in 
a CCD camera, or inferred from spectroscopy of the vibra- 
tional modes. A simpler experiment would require to set the 
the anisotropy parameter within the instability regime given 
by the displacements in Eqs. (14) and (13), which leads to 
(Kc,i — Kx)/Kcj ~ Vy ~ 10 _2 and/z//^ 10~ 2 . The linear para- 
magnet would be directly unstable towards to zigzag AF with- 
out the need of adiabatically tuning the transverse field g(t), 
and the demanding accuracy over the trap frequencies. 

Conclusions and Outlook- A sensible QS must address 
questions that are difficult to assess by other analytical or nu- 
merical methods. In this work, we have proposed a trapped- 
ion QS that fulfills this requirement. In particular, in the 
regime where non-adiabatic effects of the zigzag distortion 
become relevant, the complexity of the many-body model in 
Eqs. (3) and (6) compromises the efficiency of existing nu- 
merical methods. Besides, this QS may address the effects of 
the inhomogeneities, the long-range dipolar tail of the spin- 
spin interactions, and the dynamics across such a magnetic 
structural quantum phase transition. We emphasize that the 
incorporation of all these effects make our QS of the utmost 
interest, which may also find an application in the context of 
other Wigner crystals [36], such as electrons in quantum wires 
or liquid helium. 
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SUPPLEMENTARY MATERIAL 

1. Phonon condensation in the linear- to-zigzag Coulomb 
crystal 

In this Appendix, we present a detailed discussion of the vi- 
brational modes of an inhomogeneous Coulomb crystal close 
to the linear-to-zigzag SPT. The description for the homoge- 
neous case, which applies to ring traps or to the center of a 
linear trap in the thermodynamic limit, can be found in [1]. In 
this Appendix, we focus on the inhomogeneous crystal. 

Let us recall that the atomic plasma consists of N ions of 
charge e and mass m coupled by means of the Coulomb inter- 
action and confined by an effective quadratic potential 

v = t E f^+^ra' (16) 

At low temperatures, the ions self-assemble in a crystalline 
structure due to the balance of the trapping force and the 
Coulomb repulsion, namely V r .V = 0. The solution of this 
system of equations yields the equilibrium positions rf = l z f 9, 
where l z = (e 2 /m(D 2 ) 1 / 3 is a unit of length, and determines the 
particular geometry of the ion crystal. In this Appendix, we 
consider the regime of a linear ion crystal [Fig. 3(a)], and ex- 
pand the potential (16) around the equilibrium configuration 
Tj = l z (^e z + qi x e x + qi y e y + qi z e z ), considering oscillations of 
a small amplitude q\ a <C . The Taylor-expanded potential 
becomes V ~ Vb + Vh + V &9 where Vb stands for the configura- 
tional energy of the linear chain, Vh represents the harmonic 
approximation containing quadratic terms, and V a is the anhar- 
monic contribution due to the remaining non-linearities. The 
harmonic part can be written as follows 

Vh = J^^ma%lty¥q ia q ja (17) 

a hi 

where we have introduced the vibrational couplings 

(18) 

which depend on the anisotropy parameters K a = (co z /co a ) 2 , 
c a = 1 — 38 az , and 8 mn stands for the Kronecker delta. The 
corresponding Hamiltonian H = Yia^z /^ m )Pia + can be 
diagonalized by a canonical transformation (see e.g. [2]), 

Qua = Y^^qia, Pna = ^^Pia, (19) 

i i 

where the orthogonal matrices jtfi? fulfill 

%Ji^Jtf m = ^Km, 1 £^ i a n^"n = S i j, (20) 
ij n 

such that yield the normal-mode frequencies CO na = 
®a\j^n displayed in Fig. 3(b). In this figure, we represent 
the radial phonon modes for a chain of N = 30 ions, such 




Figure 3. Radial phonon modes for a linear ion chain: (a) 

Equilibrium positions of linear Coulomb crystal with N = 30 ions 
and (Oy = 30co z , K x < Kq. The ions self-assemble in an inhomoge- 
neous string along the z-axis, such that the distance between nearest 
neighbors is minimized at the center of the trap. Also show the laser 
wavevector 1*l within the xy-plane. (b) Vibrational frequencies co na 
for the radial modes a = x, y of a linear ion chain. As K x — >• Kq, the 
frequency of the lowest-energy mode tends to zero (0\ x — > 0. Also 
shown a scheme of a laser beam, which is detuned with respect to the 
electronic transition such that ft>2 — (Oq is close to the blue- sideband 
resonance o>2 — C0q « co ny . The large anisotropy C0y ^> co x makes the 
coupling of this laser to the vibrational modes along x negligible. 



that (Oy = 30(O z and (D x is varied. According to these normal 
modes, the harmonic Hamiltonian can be expressed as follows 

H = Y.naWV 2m ) P na + ^?GL) 5 where we have introduced 
the relevant coupling strength 

r„ =m<o„l, ^1-— ^ j. ,21, 

Note that for the so-called axial modes c z = — 2, this parameter 
is constrained to be positive r z n > 0. The same occurs for the 
radial modes along the strongly-confined y-axis. Even if c y = 
1 in this case, the vanishingly small anisotropy K y <C 1 leads to 
r^n > 0. In both cases, the additional non-linearities of V a can 
be neglected, and the vibrational Hamiltonians correspond to 
a collection of uncoupled harmonic oscillators 




where a n ,b n are the phonon annihilation operators. 

A different situation occurs for the radial modes along the 
x-axis, where c x = l, and K x < 1 is not that small. In this case, 
there might exist a certain regime where r^ s — » for a partic- 
ular vibrational mode n s G {1 • • -N}. This gapless vibrational 
mode, the so-called soft mode, is the paradigm in the Landau 
theory of SPTs [3]. As occurs in the homogeneous case [1], 
the soft mode for the inhomogeneous linear-to-zizag transi- 
tion corresponds to the zigzag mode n s = 1, which becomes 
gapless a the critical point rj s = for K x = K c [Fig. 4(a)]. In 
this situation, one is obliged to consider the anharmonic terms 

i^j iD Ki —Zj\\n m J 

(23) 

which couple the soft phonon mode to the rest of the vibra- 
tional modes nj^n s , leading thus to a complicated many-body 



7 



problem. However, the higher-energy modes can be integrated 
out, yielding a set of renormalized parameters for the soft 
mode that incorporate the effects of finite temperatures [4]. 
Since we are interested in the T = case, we neglect these 
corrections, and simply consider the effect of the quartic term 
on the soft zizag mode 

v*=\«&, < = ^££ ^# - (24) 

* * i tfi \ z i ~ z j\ 

From this expression, it is evident that u x n ^ > 0, which guaran- 
tees the stability of the theory for the soft mode in the regime 
r^ s < 0. The effective description of the SPT corresponds to 
the so-called (f) 4 scalar field theory 

/2 

#s = —Pj; +hr5(£ x + \< Qt r, (25) 

such that the broken symmetry phase rj s < 0, > yields 
the zigzag configurations of the ion crystal. Note that, as 
customary in bosonic interacting models of this type, this 
broken- symmetry phase can be understood as a consequence 
of the macroscopic occupation of the soft mode (Qn S: x) ^ 
(i.e. phonon condensation). Two comments are now in order. 
Firstly, since the model corresponds to a real scalar field, this 
phonon condensation is not associated to the phenomenon of 
superfluidity familiar to trapped gases of weakly-interacting 
bosonic atoms. Secondly, the usual long-range order associ- 
ated to the condensed state only takes place at T = 0, so this 
state of matter is rather referred to as a quasi-condensate. 

Finally, let us address the accuracy of this description by 
comparing its predictions to the estimates for inhomogeneous 
Coulomb crystals in [5], which are based on a different ap- 
proach based on the theory of prolate spheroidal plasmas. The 
critical point implicitly given by Eq. (21) is 

i K-^f 

K c = 2LL 1^0 ^0.3 (26) 
Z i rfi \ Z i ~ Z j\ 

which has been represented in Fig. 4(b), and compared to the 
analytical estimate in [5]. This critical point marks the on- 
set of the linear-to-zigzag transition, which takes place at the 
center of the trap and progressively follows to the edges. The 
agreement of both predictions justifies the validity of the re- 
sults presented in this Appendix. 



2. Inhomogeneous structural phase transition and coupled 
double-well oscillators 

In the previous Appendix, we discussed the onset of the 
linear- to-zigzag SPT for inhomogeneous ion crystals. The 
theory is based on the soft zigzag mode described by an effec- 
tive (f) 4 model, which pinpoints the origin of the instability at 
the center of the trap, but does not account for its propagation 
towards the edges. In order to describe this effect, we should 
resort to the local ion vibrations rather than to the collective 
phonon modes, but exploiting our knowledge about the zigzag 




Figure 4. Mode softening for a linear ion chain: (a) Parameters 
as a function of the anisotropy K x for a linear Coulomb crystal 
with N = 30 ions. As K x — >> K c = 1 • 10 -3 , rj s — >> for the lowest 
frequency mode, which corresponds to the zigzag distortion (see the 
mode ^i ns in inset for k x = Kq/2). (b) Scaling of the critical point 
as a function of the number of trapped ions, which follows from 
the solution of Eq. (21) (yellow dots). This parameter is compared to 
the estimate = 64(log(3Mti /2 3 / 2 ) - l)/9N 2 xi (red solid line) [5], 
where x\ = (7£(3)/2) 1//2 , and £(s) is the Riemann zeta function. 



distortion being responsible of the SPT. This can be achieved 
by the ansatz qi x = (—l) l Sqi x introduced in [6] for the ho- 
mogeneous case, which we use here for the inhomogeneous 
crystal. Note that Sqt x is a slowly varying transverse dis- 
placement, which allows a gradient expansion Sqj x w 8qi x + 
(?j -$)di+iSq ix , where djSq ix = (Sq ix - Sq jx ) <C Sq ix . Note 
that this ansatz is inspired by the actual shape of the zigzag 
mode (see inset of Fig. 4(a)), which corresponds to an al- 
ternating distortion modulated by a slowly varying envelope. 
Let us also remark that the above gradient expansion is only 
valid away from the edges, where the inhomogeneities be- 
come much larger, and higher powers of the gradient might 
be required. However, it suffices to describe how the instabil- 
ity proceeds towards the edges. The Taylor expansion of the 
Coulomb interaction gives rise to 

V x = £ (f 54 + \ifMi)) + 1 ^(djdq^) 2 + AV X , 

(27) 

where AV^ accounts for the effect of the remaining vibrational 
modes. Here, we have introduced the characteristic couplings 
of Eq. (4) in the main text, which depend upon the following 
function 

W W = L Is0_s0|n ' (28) 

where nGZ. Neglecting the effects of the remaining vibra- 
tional modes, the Hamiltonian H x — Y,i(l^/2m) pj x + V x corre- 
sponds to a set of anharmonic oscillators with site-dependent 
couplings 

rf = m(0 2 x l 2 z (1 - \k x ^{3)) , ii? = m(0 2 x l 2 z (^K x Q(5)j , 

(29) 

which clearly resemble the parameters of the zigzag (f) 4 model 
in Eqs. (21) and (24). In addition, one finds a nearest-neighbor 
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Figure 5. Local double- well structural phase transition: (a) Pa- 
rameter rf as a function of the ion position for different anisotropics 
K x in a linear Coulomb crystal with N = 30 ions. As K x — > K^n/2> 
~~ ^ ^' S1 g na l m g me onset of the SPT at the center of the trap, 
(b) Comparison of the inhomogeneous critical point Kqj given by 
Eq. (31) (yellow dots), with the analytical estimate [5] that predicts 

kc :i = Kc/(1 - A/c(z°) 2 ) 3 , with Ak = \6/(9N 2 ^ /3 x 4 / 3 ) (red line). 
We show the clear agreement between both predictions. 



coupling strength 

Kfj = matt (l^jr^J 8j»u (30) 

leading to a collection of coupled double- well oscillators. 

By numerically evaluating Eq. (30), one finds that the oscil- 
lator coupling fulfills > 0,Kfj > 0, and thus the lower en- 
ergy configuration tries to minimize the gradient. Assuming 
that (Sqj X ) = (Sqi X } in the bulk of the ion chain, the broken- 
symmetry phase (Sqi x ) ^ occurs when r^ 2 < leading to 
the double- well potential. However, the site-dependent cou- 
plings induce an inhomogeneous nucleation that originates at 
the center of the trap [Fig. 5(a)]. The inhomogeneous critical 
point is given by Kqj = 2/^(3), namely 

g -i [(-!)' -(-IV] 2 m 
A m \ z i~ z j\ 

which increases as one proceeds to the edges of the ion crystal. 
Let us remark that such an inhomogeneous critical point (yel- 
low dots) has been compared to the analytical estimate based 
on spheroidal plasmas (red solid line) [5] in Fig. 5(b). The 
clear agreement justifies the validity of the above ansatz. 

We have evaluated numerically the remaining parameters 
of Eqs. (29) and (30). In Fig. 6(a), we represent the parameter 
associated to the quartic contribution, which becomes more 
important in the vicinity of the critical point. Note that, al- 
though Sqf x <C Sqj x , the quartic contribution becomes rele- 
vant when K x — )► K c where U{ x ^> r iX . Let us also note that the 
oscillator couplings Kfj also become more relevant close to 
criticality (see Fig. 6(b)). However, they are several orders of 
magnitude smaller than the remaining parameters. This justi- 
fies our treatment of the spin-Peierls phenomena, whereby we 
have neglected their effect in a first approximation. 




Figure 6. Inhomogeneous parameters: (a) Parameter u\ as a 
function of the ion position for different anisotropics K x in a linear 
Coulomb crystal with N = 30 ions, (b) Parameter K* i+V 



3. Generalized XY model with Dipolar Interactions 

In this Appendix, we present a detailed derivation of the 
dimerized quantum Ising model (DQIM) responsible of the 
spin-Peierls transition. Let us start from the laser-ion inter- 
action in Eq.(2) of the main text, which we rewrite here for 
convenience 

H L = -^(^ / a+ + ^ / *Oi-)(e ik ^ r '- i ^ +i ^ + e - ik ^ +i ^-^). 
2 i,i 

(32) 

We now express the position of the ions in terms of the collec- 
tive phonon modes (22), namely r ; = l z £- e z + 5r ; , where 

= l z q ix e x + £ J?_- (a n + a\) + £ ~f== ip n + b\) . 

(33) 

As argued in previous sections, in order to treat the inhomo- 
geneous SPT, it is more appropriate to consider the vibra- 
tions along the x-axis in the local basis qi x . Note that the 
laser wavevectors are aligned within the xy-plane [Fig. 3(a)], 
so that the spins are not coupled to the axial phonons b ni b\, 
which shall be neglected. Besides, the large difference be- 
tween the radial trapping frequencies, co x <C C0y, implies that 
the laser beams in resonance with the vibrational sidebands 
along the j-axis will be out of resonance with respect to the 
sidebands along the x-axis [Fig. 3(b)]. Note that the reso- 
nance with higher sidebands can also be neglected by con- 
trolling the trapping frequency (Q x and the orientation of the 
laser wavevector kL, so that the coupling of the spins to the 
vibrational modes along the x-axis can be ignored. There 
is however one exception, the soft phonons can condense 
in the zigzag mode as a consequence of the SPT, affecting 
thus the laser-ion Hamiltonian. Hence, we must consider 
H o = \ Li (dotf + Li (O n ai a n + where co n = (Q ny , and 

(34) 

according to Eqs. (29) and (30). We now transform the 
laser-ion Hamiltonian to an interacting picture with respect 
to Hq = ^cr? + Y, n ^ n (^n a n- After a Taylor expansion of 
the Hamiltonian (32) for small Lamb-Dicke parameters r\\ n = 
ki y /^/2m(O n < 1, we neglect the fast-oscillating terms from 
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the resulting spin-phonon Hamiltonian by using a rotating- 
wave approximation (RWA). This RWA holds when the laser 
beams are tuned to the red (0\ ~ cOq co v , and blue o>2 ~ C0o + 
0)y sidebands, and fulfill the following constraints £2/ <C COo, 

T7i^i<C|0)i-(0)o + ^)|, T] 2n ^ 2 <|6>2-(^-^)|. (35) 

In this regime, the laser-ion Hamiltonian becomes a combina- 
tion of the red- and blue-sideband spin-phonon couplings 

(36) 

where we have introduced the sideband coupling strengths 
&\ n = \£l x r\i n J{ in ^\^ n = ^D. 2 ri2n^inQ i<h , and the detun- 
ings 8 m = CO] - ((Oq- CD n ),8 nb = o>2 - ((Qo + co n ). Let us re- 
mark that in this spin-phonon Hamiltonian, the zigzag dis- 
placement that pinpoints the structural phase transition of the 
ion crystal is also considered in qi x = ( — l) l 8qi x , such that 
the commensurability parameters are 6 r — k\ x l z ,Q\> = k2 X l z . 
By making the sideband detunings opposite to each other 
8m = —8 n b ='• 8n, it is possible to move to a picture where 
the phonons absorb the time-dependence, and the laser-ion 
Hamiltonian becomes 

Hi = £ 8 n ala n + £ (^ n a n + H.c.) , (37) 

n in 

where we have introduced the following spin operator 

fin = (^)*a,-e- i0 ^ +^a+e +i ^^. (38) 

The current formulation of this generalized spin-phonon 
Hamiltonian (37) is amenable of performing perturbation the- 
ory in the regime <C 8 n . Note that in this limit, the 
phonons along the v-axis only get virtually excited and it is 
possible to integrate them out to obtain an effective Hamilto- 
nian that only involves the spin degrees of freedom. Such type 
of perturbative expansion can be performed systematically by 
a polaron-type transformation 




At this point, we perform a mean-field approximation for 
the zigzag mode, so that qi x = (qi x ) can be treated as a c- 
number. The transformed Hamiltonian can be obtained to any 
order of perturbation theory by using the identity e s H^e~ s = 
H[ + [S,Hl] + ± [S, [S,HL]] + • • • . To lowest order, we find 
UH{U^ « Y,n 5n<Ai a n +#eff> where 

ff eff =11 T^jn + L ^nmOf , (40) 
ij n n inm 

and we have introduced the following phonon operator 

i /■ <^^IYI ^^\TYl / 2 2 \ / "f" "f* \ 

(41) 



From this last expression, it is clear that by setting the laser 
parameters according to the following constraint 

a\k\ y = a 2 2 k 2 2r (42) 

the residual spin-phonon operator vanishes, and we obtain 
an effective model of interacting spins. In particular, this 
constraint is satisfied for a pair of counter-propagating laser 
beams ki = — k 2 =: k with identical Rabi frequencies £l\ = 
Q.2 =' By substituting the generalized spin operators 
in (38), we obtain the following effective spin Hamiltonian 

¥J v ' ' ' 7 

(43) 

where the coupling strengths are those written in Eq. (8) of the 
main text. This Hamiltonian corresponds to a generalization 
of the famous XY model [7], where in addition to the well- 
studied anisotropy between the Jf* and Jjj couplings, there 
are additional crossed couplings JfJ^Jj which have remained 
largely unexplored. These couplings are expressed in terms 
of the relative phase between the lasers 0_ = <j>\ — 02, and the 
strength 

Q2 

Ji J = ~ L JZt rinTlm^in^jm • (44) 

Note that in the limit of tight transverse confinement K y = 
( (O z /(Oy) 2 <C 1, the above couplings can be expressed as 

•Jeff ^L^v^y / 8y\ 

J -'=wf^f J " = Wv + %)'"'• <45) 

where we have introduced the bare detuning Sy = (0\ — (coo — 
CQy) = —0>2-\-((Do + (Qy). Therefore, the effective spin model is 
precisely the generalized XY model with dipolar interactions 
expressed in Eqs. (8)-(9) of the main text. 

In order to obtain the final expression of the Hamiltonian, 
note that the laser wavevectors lie in the xy-plane and are 
slightly tilted with respect to the v-axis [Fig. 1(c) main text]. 
Accordingly, we Taylor expand the spin couplings for the 
small parameter ^ = Q{8qt x ) <C 1. By considering that the 
laser beams are phase locked 0_ = 0, we obtain the announced 
dimerized QIM 

J% = 2J ij (-l) i+1 =J$. (46) 

4. Dimerized Quantum Ising Model 

a. Analytical Solution 

In this part of the Appendix, we present the analytical solu- 
tion for the energy spectrum of the DQIM in Eq. (11) of the 
main text. By restricting to nearest neighbors in a homoge- 
neous crystal, the Hamiltonian becomes 

ffdf = I/(ofof + i-(-l) < «ofof +1 + (-l) i 4ofof + i-«of), 

i=l 

(47) 
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where J > stands for the antiferromagnetic spin coupling, 
and t; , g > represent the dimerization and transverse field 
expressed in units of the spin coupling. In contrast to previ- 
ous studies [8, 9], the dimerization term in Eq. (47) does not 
commute with the bare spin couplings and thus introduces ad- 
ditional quantum fluctuations. Motivated by the derivation of 
the Hamiltonian for a system of trapped ions presented in the 
previous section, we shall focus on the limit of small dimer- 
ization £ <C 1 , where the magnetic phases coincide with those 
of the usual QIM, namely, the paramagnetic (P) and the anti- 
ferromagnetic (AF) phases. As shown below, the dimerization 
is responsible of (a) lowering the ground- state energy, and (b) 
displacing the critical point. This two-fold effect will be ex- 
ploited to show that it is possible to observe a spin-Peierls 
quantum phase transition with trapped ions. 

We diagonalize the Hamiltonian (47) via the so-called 
Jordan- Wigner (JW) transformation [10]. We note that the 
effect of alternating dimerized couplings in the XY model 
has been addressed in [9]. To the best of our knowledge, 
the dimerized crossed couplings in Eq. (47) have not been 
studied so far. The JW mapping expresses the spin operators 
in terms of fermionic operators. The spin up (down) states 
are associated to the presence (absence) of a spinless fermion 
It)/** c Il vac )z> I \)i ** l vac )z> by a non-local transformation 

g+ = c]j n ^< iC "i c i, or = c- iK ^< iC J c J Ci , of = 2c] Ci - 1, 

(48) 

where we have introduced of = (of ±iof )/2. The fermion- 
ized Hamiltonian can be expressed as follows 

#eff = L / (( 1 + i (-l)^)cJ +1 c- -2gc] Ci -cj +1 cj) +H.C, 

(49) 

which corresponds to a tight-binding model of spinless 
fermions with dimerized complex tunnelings t — J (I =bi2§), 
on-site energies £ = —2Jg, and a pairing term between nearest 
neighbors that can be interpreted as the injection of Cooper 
pairs from a neighboring superconductor with gap A = —J. 
This Hamiltonian is diagonalized in momentum space Cj = 
Y<qZ lq ^ c qi where q G [— 7T, n) stands for the crystal momen- 
tum, and we have set the lattice spacing to unity. By intro- 
ducing the Nambu spinor ^(q) = (c^,cl^,c^_ 7r ,c^ +7r ) / in a 
reduced Brillouin zone, the Hamiltonian becomes 

H= £ V(q)lU q y(q) : RBZ=[0,7T/2), (50) 
qeRBZ 

where is expressed in terms of Pauli matrices 

M q =2J (cos qT z (g>T z + sinqT z (g)Ty — 2{; cosqZy®l2 — gl2®T z ). 

(51) 

The diagonalization of this matrix yields four eigenvalues 

U q M q Ul = dmg{-£-(q),-£+(q),+£-(q),+£+(q)}, where 



£±(q) = 2jJ (g ± v / l+4£ 2 cos4) 2 + sinV 



(52) 











c 
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Figure 7. Scaling of the groundstate energy (a) Groundstate en- 
ergy of the DQIM as a function of the dimerization parameter ^ . The 
exact expression £g W (<^) (53) is represented by yellow triangles, and 
compared to the scaling £| pp (^) (55) in blue solid lines In the in- 
set, we display the error due to the approximation for ^ <C 1, namely 
|AEg| = |£g W — Eg PP \. (b) The function <p(g), which determines the 
quadratic scaling of the groundstate energy (55), is a positive-definite 
and monotonically-decreasing function. 



of the system corresponds to all the negative-energy levels 
filled |g) = n q 7q i Oq" 2 1 vac )' where |vac) stands for the vac- 
uum. In the thermodynamic limit, the ground- state energy 
£g W = — Lg G RBz( e - (tf) + e + (#)) can he expressed as follows 

IN r n I , 
= ~— / <W(S- V / T+4| I cos, ? )2 + sinV 

& Tt JO v 

(53) 

To check the consistency of this result, note that this en- 
ergy coincides with the ground- state energy of the standard 
QIM [1 1] in the limit of vanishing dimerization £ = 0, namely 
£^ IM (0) = -ljN(i+ g )g>(Q g ) m Here, we have introduced the 
complete elliptic integral of the second kind 

K 

£{0 g ) = Jj da(l-0 2 sin 2 a) 1/2 , (54) 



The associated unitary matrices define a generalized Bogoli- 
ubov transformation (y q ^ , y^ 2 , 7^,3, YqA Y = U q ^!{q) that di- 
agonalizes the fermionized Hamiltonian. The groundstate 



where 6 g = y/4g/(l -f g) 2 , and have used its properties listed 
in [12]. With these expressions, we can address the effects of 
the dimerization announced at the beginning of this Appendix. 

(a) Lowering of the ground-state energy - Considering the 
limit § <C 1 and the ground-state energy in Eq. (53), we find 

E™($) -E™(0) = \l d£'^ = - 2 -f<p(g)e, (55) 
where we have introduced the following function 

rK 

<p(g)= / dq(cos 2 q-gcosq)(l+g 2 -2gcosq)- l/2 . (56) 
./() 

By using the complete elliptic integral of the first kind 

/*! 

Jf(0 g ) = J da(l-e g 2 sin 2 ar 1/2 , (57) 
this function can be expressed in a compact form 

<Pte)= 3?( <f( ^ )(1+g)(2 ^" 1)+ ^ (0g)(1 "^ )(2g2 + 1) )' 

(58) 
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Figure 8. Displacement of the critical point: Magnetic suscepti- 
bility Xm as a function of the transverse-field ratio g. The divergence 
of the susceptibility pinpoints the critical point g c where the quan- 
tum phase transition from a paramagnet to an antiferromagnet takes 
place. The arrows indicate the flow of the critical point g c — > g c (£) 
for different dimerization parameters e {0, 0. 1 , 0.2, 0.3, 0.4, 0.5} 



The above equation (55) is precisely the expression (12) used 
in the main text of this article. To test its validity, we compare 
it to the exact expression E™(1;) in (53) for g — 0, where 
<p(0) = 7T/2 (see Fig. 7(a)). As shown in the inset of this 
figure, the agreement between both expressions is extremely 
good even for considerable dimerization parameters £ < 0.4. 
Since this regime contains the limit that can be attained with 
the ion trap QS, we can use directly the groundstate energy 
Eg W (^). From Fig. 7(b), and from the properties of the el- 
liptic integrals, <p(g) is positive-definite function that gives 
rise to the aforementioned lowering of the groundstate energy 
E™(Z)-E™(0) « -^r\(p(g)\Z 2 < 0. This property be- 
comes crucial for the spin-Peierls transition, since it will com- 
pensate the energy increase due to the structural ion change. 

(b) Displacement of the quantum critical point - As men- 
tioned at the beginning of this Appendix, the DQIM for £ <C 1 
encompasses two magnetic phases that coincide with those of 
the usual QIM, namely, the paramagnetic and antiferromag- 
netic phases. The antiferromagnetic (AF) phase has a degen- 
erate groundstate that corresponds to the two antiparallel Neel 

configurations |AF) £ {| H 1 — > , | — I h)} in the 

limit of g <C 1 , where we have introduced the up/down spins 
in the x-basis |±) = (| t) =M i))/V2. This degeneracy is re- 
lated to the invariance of the Hamiltonian under a global spin 
inversion U% 2 = fl of > which also applies to the DQIM. Con- 
versely, for g ^> 1, the paramagnetic (P) phase has a unique 
groundstate where all the spins are parallel to the transverse 
field |P) = | tt " ' t)- As the transverse field is decreased, the 
paramagnetic groundstate evolves towards one of the possible 
antiferromagnetic states, breaking spontaneously the symme- 
try at the quantum critical point g = g c . At this symmetry- 
breaking point, the energy gap between the groundstate and 
the lowest-lying excitations vanishes. 

For the QIM, the critical point is located at g c — 1, as im- 
posed by the self-duality of the Hamiltonian under the order- 
disorder transformation [13], /if = Ylj<i a h Mf = a f-\ a f^ 

H({o},g) = -gH({Li},g- 1 ) = -ZMtiri+i -*"Vf), 

(59) 

which implies that the energy gap A(g c ) = = A(l/g c ) only 



vanishes at g c = (g c )~ l , such that the critical point gets 
locked to g c = 1. Note that this self-dual symmetry is ex- 
plicitly broken by the dimerization (47), which lifts the con- 
straint and allows the quantum critical point to attain differ- 
ent values. In Fig. 8, we represent the magnetic susceptibility 
Xm * Xm = ~d 2 Eg/ dg 2 8LS3. function of g for different dimer- 
izations & e {0,0.1,0.2,0.3,0.4,0.5}. The divergence of the 
susceptibility marks the critical point of the QPT. From this 
figure, it is evident that the critical point g c = 1 flows towards 
higher values g c (%) as the dimerization increases. 

In order to find analytically the location of the displaced 
critical point, we shall use the explicit expression for the en- 
ergy bands (52). The lowest-lying excitations occur at q = 0, 
such that their energy gap scales as A oc \g— ^\ +4^ 2 |. The 
divergence of the correlation length typical of critical phe- 
nomena, which implies A(g c ) = 0, leads to 

g c ^| c (0 = V / T + 4F, (60) 

which is precisely the result in Eq. (13) of the main text. 
Hence, the critical point is displaced towards higher values 
of the transverse field. This effect will also turn out to be es- 
sential for the quantum nature of the spin-Peierls transition. 

b. Numerical Analysis 

The increased complexity of the model brought by the 
dimerization demands a careful assessment of the validity of 
the above results. In this part of the Appendix, we make a 
direct comparison between our analytical derivations and the 
numerical results obtained by the density-matrix renormal- 
ization group (DMRG) method [14]. DMRG algorithms are 
considered nowadays as the most efficient approach to one- 
dimensional QMBS, and therefore ideally suited to explore 
the strongly-correlated spin model (47). We focus on the 
groundstate properties of the model, and provide numerical 
evidence that clearly supports the validity of the expressions 
for the energy (53), and its scaling with the dimerization (55). 

We use the finite-system DMRG algorithm for open bound- 
ary conditions. In order to test the accuracy of our algo- 
rithm, we have compared the obtained groundstate energies 
£E>mrg tQ the resu i ts f a Lanczos algorithm ^Lanczos for 

the exact diagonalization of small spin chains of length L £ 
{10, 12, 14, 16, 18,20,22,24}. Note that in order to achieve 
^max = 24, we have optimized the Lanczos algorithm using 
the bipartite nature of the spin model (47). By dividing the 
chain in two blocks, it is possible to lower the computational 
cost of matrix-vector operations, thus optimizing the perfor- 
mance of the Lanczos algorithm. For the DMRG algorithm, 
it suffices to keep m = 20 states of the reduced block density 
matrices, and considered N s = 1 full sweeps of the renormal- 
ization procedure. In Fig. 9(a), we represent both energies for 
the DQIM (47) after setting § = 0.1, g = 0. The good agree- 
ment serves as a testbed for our DMRG algorithm, which we 
now use to test the analytical prediction (53) valid for the ther- 
modynamical limit L — >• oo. 

Let us remark that the exact diagonalization cannot be car- 
ried far beyond L max = 24. In clear contrast, the DMRG al- 
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Figure 9. DMRG calculation of the grounds tate energy: (a) Com- 
parison of the groundstate energy of the DQIM obtained by exact 
diagonalization methods based on a Lanczos algorithm Eg 3110208 (£) 

(yellow dots), and the DMRG method Ef MRG (%) (dashed grey line). 
We analyze § = 0.1, g = for small chains with L < 24 sites. In 
the inset we represent the error \AE f g \ = |^Lanczos _^DMRG| (blue 
crosses), (b) Comparison of the analytical estimate for the ground- 
state energy E 3 g w (^) (blue solid line) in the limit L — » °o, and the 
DMRG result £^ MRG (<!;) (red diamonds) for longer chains of up to 
^max = 400 sites, and %=0A,g = 0. In the inset we represent the 

(green crosses). 
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5. Spin-Peierls Quantum Phase Transition 

In this Appendix, we discuss the possibility of adjusting 
the parameters of the spin-phonon model, such that the mag- 
netic and structural phase transitions occur simultaneously as 
a consequence of the quantum fluctuations introduced by the 
transverse field g. Let us recall that under the adiabatic and 
homogeneous approximations for the (f) 4 model, the lowering 
of the groundstate energy (55) can be rewritten as follows 

E g ^)-E g (0) = - 2 ^ 1 L^l (61) 



Since this expression is quadratic in the zigzag displacements, 
we must modify the parameters of the scalar field theory (29) 
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In the homogeneous approximation, whereby the coupling be- 
tween the neighboring double- welled oscillators is neglected, 
the SPT occurs when rf < 0, u\ > 0. Hence, it becomes evi- 
dent that the critical point is displaced towards a smaller value 
of the ratio between the trapping frequencies 
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Figure 10. DMRG energy scaling law: Comparison of the ground- 
state energy obtained by the DMRG method ^ MRG (^) (red dia- 
monds) with the Jordan-Wigner solution E™(t;) (yellow triangles) 
and the approximate expression £| pp (^) (blue solid line). In the in- 
set we represent the corresponding errors \AE r g \ = |£| pp — i^ MRG 
(blue-red crosses), and \ AE'\ = |£g W — £J? MRG | (yellow-red crosses) 



Hence, the lowering of the energy due to the magnetic order 
compensates the energy increase due to the zigzag distortion, 
and the linear chain becomes unstable towards the distorted 
lattice for %^ < k x < K^-. 

The idea now is to tailor the parameters of the model in or- 
der to make the critical points of the magnetic (60) and struc- 
tural (63) phase transitions coincide. This occurs for 



gorithm can be applied to longer chains, allowing us to con- 
sider L e {40,80,120,160,200,240,280,320,360,400}. In 
Fig. 9(b), we represent the obtained energy £^P MRG , and 
compare it to the Jordan-Wigner estimate £"g W (53) for § = 
0.1, g = 0. In this case we have set m = 50 and N s = 2. Note 
that the analytical prediction gets more accurate as the ther- 
modynamical limit is approached (see the inset of Fig. 9(b)). 



2J 



6(3) 



71 



o 2 <p(gcY 



(64) 



which is equivalent to the condition for the axial trapping fre- 
quency in the main text (15). By substituting in Eq. (62), we 
find that 



ff=^/fSi(3)(l-J^). 



(65) 



Finally, we have checked the scaling of the groundstate en- 
ergy with the dimerization parameter (55). We consider a spin 
chain with L = 200 sites, and use the finite- size DMRG al- 
gorithm for m = 50 and N s = 2. In this case, we vary the 
dimerization ^ G [0,0.5] for g = 0, and study the behavior of 
the groundstate energy (see Fig. 10). The agreement shown 
in this figure, even for relatively large dimerizations, provides 
sufficient evidence to support the claims of this work. 



Due to the monotonically-decreasing nature of the function 
<P(gi) < <P(g2), for gi > g 2 , (see Fig. 7(b)), we find that 
g > g c rf > 0, leading to a disordered linear paramagnet, 
whereas g < g c rf < yields the ordered zigzag antifer- 
romagnet. Hence the spin-Peierls transition is triggered by 
quantum fluctuations alone (i.e. spin-Peierls quantum phase 
transition), as outlined in the scheme of Fig. 1 in the main 
text. 
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6. Spin-Peierls Quantum Simulator in Perspective 

In order to assess the scope and utility of the proposed quan- 
tum simulator, we discuss in detail some of the experimental 
and theoretical efforts to understand the spin-Peierls instabil- 
ity in the context of quantum many-body physics. 

Experimental perspective - The first experimental evidence 
of the spin-Peierls instability was found in the organic molec- 
ular crystal TTF-CuS4C4(CF3)4. There, unpaired spins s = 
1/2 localized at neighboring donor ions TTF + interact via 
super-exchange through the acceptor CuS4C4(CF3)4, and give 
rise to an alternating antiferromagnetic Heisenberg model that 
depends on the lattice dimerization. The opening of the spin- 
Peierls gap was inferred from the sharp isotropic decrease of 
the magnetic susceptibility below T c « 12K [15], and from an 
anomaly in the magnetic contribution to the specific heat [16]. 
Besides, the associated dimerization of the lattice was directly 
inferred from X-ray measurements of new Bragg peaks [17], 
the intensity of which gave evidence for the condensation of 
a soft phonon mode. Finally, the clear-cut confirmation of 
the magnetic origin of the structural phase transition was ob- 
tained by applying an additional magnetic field. In the low- 
field regime, neutron scattering experiments displayed a shift 
of the critical temperature [18]. Additionally, magnetization 
measurements showed that a new phase transition takes place 
in the high-field regime [19, 20], and leads to an incommen- 
surate lattice distortion. Let us finally note that spin-Peierls 
transitions were also observed for other organic materials, 
such as TTF-AuS 4 C4(CF 3 )4 (T c « 2. IK) [21], and MEM- 
(TCNQ)2, which displayed both an electronic Peierls transi- 
tion (T c ^335K) and a spin-Peierls instability (T c «18K) [22]. 

The discovery of a spin-Peierls phase transition in CuGeC>3 
(T c ^14K) [23], an inorganic compound, paved the way to 
new experiments due to the availability of larger crystals with 
a higher quality. In this case, the spins s = 1 /2 of localized 
electrons in the Cu 2+ ions interact via super-exchange through 
the O 2- ions, and also lead to an antiferromagnetic Heisen- 
berg model that depends on the lattice dimerization. Hence, 
similar effects were observed in the experiments, such as the 
drop of the magnetic susceptibility [23], dependence on ad- 
ditional magnetic fields [24], or X-ray and neutron scattering 
measurements of the secondary Bragg peaks [25]. In addi- 
tion, this new material allowed for new possibilities, such as 
the measurement of the spin-Peierls gap and the low-energy 
dispersion relations via inelastic neutron scattering [26], or 
the accurate estimation of the modified lattice constants by 
neutron diffraction [27]. Let us remark that no evidence of 
the soft-phonon condensation was found in these experiments. 
Another difference with respect to the organic compounds is 
the presence of magnetic frustration in the underlying Heisen- 
berg model, which is based on the analysis of magnetic sus- 
ceptibility [28] and magnetostriction experiments [29]. 

Let us finally note that other inorganic materials, such as 
the transition metal oxide TiOCl [30], have also shown clear 
evidence of an underlying spin-Peierls transition. This mate- 
rial, which was originally proposed to realized a resonating 
valence bond state, has turned out to be another example of 
an antiferromagnetic Heisenberg magnet that displays a spin- 



Peierls transition at a much higher temperature T c « 67 K. 

Having described these seminal experiments, it is now pos- 
sible to understand the interest of the proposed quantum sim- 
ulator from a experimental perspective. First, we note that all 
of the above spin-Peierls instabilities occur at a finite tempera- 
ture T c > 0. To the best of our knowledge, the observation of a 
spin-Peierls transition only driven by quantum fluctuations is 
still an open problem, which could be addressed with the pro- 
posed quantum simulator. Second, we also note that the above 
transitions are well described by an antiferromagnetic Heisen- 
berg model. Hence, another open problem is the observation 
of spin-Peierls phenomena for materials described by other 
microscopic models. Unfortunately, the magnetic interaction 
in the above materials is determined by a super-exchange that 
cannot be experimentally controlled. In contrast, our quantum 
simulator would allow for the possibility of tailoring the mag- 
netic interactions experimentally. Although we have consid- 
ered an Ising-type coupling, we emphasize that different mod- 
els could also be explored by exploiting other laser-induced 
spin-phonon couplings. Third, an additional property of the 
trapped-ion quantum simulator is the capability of perform- 
ing very accurate measurements at the single-particle level. In 
comparison to the global measurements in the above experi- 
ments, trapped-ion technologies would allow for direct mea- 
surements of the local magnetization, or the distance depen- 
dence of two-spin correlators. Besides, the crystal structure, 
and thus the onset of the lattice dimerization, can be directly 
imaged on CCD cameras. Finally, the possibility of control- 
ling the number of ions would allow for a very interesting 
transition from the few- to the many -body scenario. 

Theoretical perspective - The first theoretical work on the 
spin-Peierls instability [31] conjectured the onset of a struc- 
tural phase transition in molecular crystals, which would be 
caused by collective magnetic interactions. This was followed 
by the first attempts towards a microscopic theory, both for 
the alternating Heisenberg [32, 33] and XY [34] magnets. 
In the latter, the spin model can be mapped onto free spin- 
less fermions, and thus connected with previous work on the 
fermionic Peierls transition (see e.g. [35]). The alternating 
Heisenberg model [33] was initially treated within mean-field 
theory (i.e. Hartree-Fock approximation), such that the results 
of the XY model could be used after a simple renormalization 
of the spin-spin interaction strength. In particular, this mean- 
field theory predicts the sharp decrease of the magnetic sus- 
ceptibility below T c , and was used satisfactorily to explain the 
experimental results on the organic crystals described above. 
The surprising accuracy of the mean-field predictions could 
only be understood in the light of the Luttinger-liquid treat- 
ment [36], which predicted only small corrections for the mea- 
sured macroscopic quantities, with the possible exception of 
the experiments with strong magnetic fields [20] . 

The challenge to theoretical theories beyond the mean-field 
treatment changed considerably with the discovery of the in- 
organic spin-Peierls sample CuGeC>3. First, the underlying 
microscopic model must incorporate a certain amount of frus- 
tration caused by next-to-nearest neighbor interactions [28]. 
We note that the interplay of frustration and quantum fluctua- 
tions is responsible for a complex phase diagram already for 
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the bare spin model [37]. In the complete spin-phonon sys- 
tem, the spin-Peierls instability and the frustration constitute 
two different mechanisms for opening an energy gap, which 
lead to two different types of phase transitions (i.e. second or- 
der and Kosterlitz-Thouless, respectively). Second, the lack of 
experimental evidence supporting a soft-phonon mode ques- 
tions the validity of the so-called adiabatic approximation for 
the crystal vibrations. In most of the theoretical approaches 
above, the lattice dynamics is neglected by only considering 
a static elastic distortion. Some of the first attempts to ac- 
count for the lattice dynamics considered the opposite diabatic 
limit [38], where the phonons, which are much faster than 
the spin excitations in this regime, contribute to a phonon- 
mediated interaction. In fact, this mechanism can be also re- 
sponsible for the long-range frustrating interactions found in 
CuGeC>3. Let us note, however, that the spin-spin coupling 
in this material has roughly the same order of magnitude as 
the phonon frequency, which seems to point out to a regime 
which is neither adiabatic nor diabatic. In this regime, one 
must resort to numerical tools, such as optimized exact diag- 
onalization [39], density matrix renormalization group algo- 
rithms [40], or quantum Monte Carlo [41]. These numeri- 
cal results showed that the non-adiabaticity leads to important 
renormalization of the spin-Peierls gap and the lattice dimer- 
ization with respect to the static theories. More importantly, it 
can also destabilize the spin-Peierls phase for strong-enough 
spin-phonon couplings via a Kosterlitz-Thouless transition. 

After this description, we can now discuss the interest of 
the proposed quantum simulator from a theoretical point of 



view. First, the trapped-ion setup is an excellent playground 
where to study strong deviations from mean-field theories. 
Second, the quantum simulator incorporates naturally the ef- 
fects of magnetic frustration, since the spin interactions de- 
velop a dipolar long-range. We note that the degree of frus- 
tration, namely the ratio between nearest and next-to-nearest 
couplings, could be tailored by slightly tilting the laser beams 
in such a way that the effective wavevector also has a compo- 
nent along the axis of the chain. Hence, the interplay of quan- 
tum fluctuations and frustration can de actively designed in 
the quantum simulator. Finally, and most notably, the degree 
of adiabaticity can also be controlled. We have focused on a 
setting where the trapping frequencies lie close to the critical 
point of the structural transition. In this regime, the frequency 
of the soft-phonon mode lies below the effective spin-coupling 
strengths, and the adiabatic approximation holds. However, 
the quantum simulator could in principle also explore the di- 
abatic regime and the intermediate one by simply changing 
the trap frequencies. From our point of view, this is a very 
important property of the quantum simulator, since analyti- 
cal approaches do not exist in the intermediate regime. Be- 
sides, numerical methods may even lead to some controversy 
when the dispersion of the phonon modes is taken into ac- 
count [42]. We would like to note that the efficiency of the 
numerical methods is likely to get degraded when the (j) 4 non- 
linearities are considered. A priori, the proposed spin-Peierls 
quantum simulator has the potential of exceeding the power 
of numerical methods on classical computers, which should 
be the final goal of any sensible quantum simulator. 
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